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Abstract 

^ , We present two pairs of Y(s/2) Yangian symmetries for the trigonometric 

and hyperbohc versions of the Hubbard model with non-nearest-neighbour 
! hopping. In both cases the Yangians are mutually commuting, hence can be 

o 
o 



combined into a Y(sZ2)©Y(sZ2) Yangian. Their mutual commutativity is of 
dynamical origin. The known Yangians of the Haldane-Shastry spin chain 
^ ' and the nearest neighbour Hubbard model are contained as limiting cases of 



our new representations. 



Yangian quantum groups were introduced by Drinfel'd more then ten years ago. 
His original intention was to put the algebraic structure underlying the Yang-Baxter 
exchange relations with rational i?-matrix, responsible for the integrability of the most 
prominent integrable systems, into the mathematically more conventional context of 
Hopf algebras. Yangian quantum groups and their representation theory are thus inti- 
mately connected with the classification of integrable quantum systems. 

*e-niail: frank@monet.phys.s.u-tokyo. ac.jp 
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However, during the past few years it became apparent that Yangians can also play 
a physically very interesting role as additional symmetries of integrable systems, and 
moreover, that Yangians are part of the symmetry algebra of such well studied integrable 
systems as the nearest neighbour Heisenberg model or the nearest neighbour Hubbard 
model 0, if considered on infinite lattices. These symmetries have been overlooked for 
many years, since it was unusual to deal directly with infinite systems. Instead all 
conventional approaches to integrable systems like the various Bethe Ansatze start from 
finite systems, usually under periodic boundary conditions, and the thermodynamic limit 
is only performed at a later stage of calculation. The Yangian symmetries of the nearest 
neighbour Hubbard and Heisenberg models are incompatible with periodic boundary 
conditions. For this reason they do not combine with Bethe Ansatz methods. 

Now there are two recent developments that make a utilisation of Yangians for inte- 
grable systems feasible. First, methods have been developed to deal directly with infinite 
systems. Within this so-called symmetry based approach it became possible to calculate 
for instance higher order spin correlators for the XXZ-chain [Q. Quantum groups play 
an essential role here. Second, interesting integrable systems with a finite number of 
degrees of freedom have been discovered which exhibit Yangian symmetry compatible 
with periodic boundary conditions, most prominent among these the Haldane-Shastry 
spin chain ||. Until the discovery of its Yangian symmetry the high degeneracy of 
its spectrum remained a puzzle, which is now resolved by Yangian representation theory 

i- 

Below we present a new pair of mutually commuting representations of the Y(s/2) 
Yangian in terms of Fermi operators, which form a Yangian symmetry of the Hubbard 
model with non- nearest-neighbour hopping of Gebhard and Ruckenstein The 
model itself contains the usual Hubbard model as well as the Haldane-Shastry spin 
chain as certain limiting cases, and so do the generators of its Yangian symmetry. 

The model describes itinerant electrons of spin o created by c^^ at site j of a one 
dimensional lattice. The probability amplitude for hopping between sites j, k will be 
denoted by tj^. Two electrons of different spin encountering each other on the same 
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lattice site feel an on-site repulsion U > 0. 

H = T.tjk4^<^ka + 2UY.{cl^c,^-l){c+^c,i-l) . (1) 

j,k j 

The constants here have been chosen for later convenience, and tjj := by definition. 
Throughout this letter we are using sum convention with respect to Greek indices. We 
will consider two different choices of translational invariant hopping amplitudes, tjk = 
tj-k- For 7^ let 

tn := -ish{K)sh~^{Kn) , n^O . (2) 

Then our choices are given hj k = in/N for a finite lattice of sites, and k, > for an 
infinite lattice. The energy scale has been chosen such as to give hopping amplitudes of 
absolute value 1 between neighbouring sites. The summation indices run from to A^ — 1 
in the trigonometric case, and over all integers in the hyperbolic case. The thermody- 
namic limit of the trigonometric model and the limit /t — > of the hyperbolic model 
coincide. In both cases tjk turns into — i/ (j — k). This is the true 1/r Hubbard modelQ. 
In the limit k oo, the hyperbolic model turns, up to a canonical transformation, 
which is described below, into the nearest neighbour Hubbard model. 

In order to understand the physical meaning of the above kind of hopping amplitudes, 
one has to consider the dispersion relation of the free model {U = 0) |T0|. In the 
trigonometric case we obtain 

iV /7r\ 

e(p)=Et„e'^" = -sin - (vr-p) , (3) 

n=l \NJ 

where p = 27i{m + 1/2) /N, m = 0, . . . , N — 1. In the thermodynamic limit this yields 
e{p) ^ 71 — p. The dispersion relation is linear in the first Brillouin zone, the model 
is chiral. It contains only left moving particles. The physically most interesting point 
about this chiral model is the appearance of a Mott transition at finite f/ > In 
*It may be interesting to notice, that the trigonometric and hyperbohc hopping amphtudes 
above can be interpreted as g-deformed l/r-hopping. The notion of q-deformation is defined by 
rq :— (g'' — q^^)/{q — q^^). Setting q equal to e'^ the hopping amphtudes become tjk ~ ^i/O ^ ^)q- 
The trigonometric case corresponds to q being the 7V-th root of unity, the hyperbohc case to g > 1. 
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the hyperbolic case the dispersion relation is 

-(p) = i:^^e-=2sh(.)f:^ . (4) 

n=l n=l ^l^yl^n) 

The last expression is easily recognised as being, up to a redefinition of scales, the 
logarithmic derivative of the Jacobi theta function ?94. As a function of n it interpolates 
between the sinusoidal dispersion relation of the of the nearest neighbour model and the 
saw tooth shaped dispersion relation of the 1/r model. 

The local U(l) transformation cj^j —>■ e^'^^Cj^, (pj real, does not alter the canonical an- 
ticommutation relations between the Fermi operators. The local electron densities c'j^Cj^ 
are invariant under this transformation, hence the interaction part of the Hamiltonian 
(H) is as well. This means that we can always use a U(l) transformation to modify the 
hopping term to our convenience. The modified model will be completely equivalent to 
the original one. Consider the case cpj = jir. This transformation introduces a factor of 
(— 1)-'"'^ into the expression for the hopping amplitudes and shifts the dispersion rela- 
tions by a half period. Using this transformation our conventions meet the conventions 
of Gebhard and Ruckenstein . To recover the nearest neighbour Hubbard model in its 
familiar form, we do not only have to consider k oo, but in addition the above trans- 
formation with (fj = jTT/2. This transformation removes the factor of "i" in front of the 
hopping amplitude, changes the hopping amplitude to an even function, and shifts the 
dispersion relation by a quarter period. Hence the quadratic bottom of the sinusoidal 
band is shifted to p = 0. 

There is yet another important canonical transformation, namely 

Cji , CjT ^ c+ , -U . (5) 

This transformation leaves every Hamiltonian of the form (1) with antisymmetric hop- 
ping matrix invariant, but the global spin operators and, in our case, the Yangian gen- 
erators (see below) are not. It is responsible for the doubling of the Yangian. 

The natural language for writing down the s/2 generators of the rotational symmetry 
of the Hamiltonian (|^) is, of course, in terms of spin operators, which are linear combi- 
nations of products of one creation and one annihilation operator at the same site. For 
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the formulation of our Yangian generators below it turns out to be useful to extend this 
concept to spin-like operators with indices corresponding to different sites. We arrange 
the pair of operators Cj'^Ckr in a 2 x 2-matrix labeled by spin indices a, r in the usual 
tensor product convention, {Sjk)^ '■= Cj^Ckn and then set 

S'jk '■= tr(a"S'jfc) , Sji^ := ti^Sjk) , 3°" := Sfj , Sj := Sj^ , (6) 

where the a" are the Pauli matrices, and the bar denotes complex conjugation. Our 
definition implies (5'";.)"'" = 5*^^-, {Sjf.)'^ = S^j. ^S" and 5*° are the spin density and 
electron density operators, respectively. The algebra of the operators 5*°^ is rather rich. 
We obtain the commutators 

Shr] = ^klSf^ — SmjS^i^ , (8) 

[S-k,sL] = 6^^{6MS^^-6^,Si)+ie^^^{6kiSl, + 6^,Sl) . (9) 

However, there are other relations. For the construction of our Yangian generator and 
the verification of the Yangian Serre relations below, we further need the following. 



'-'jk'^lm "T '-'jk'-'lm "r '^'-'jm'^lk 


— 4:^klSjjyi + 2SlmSjk , 


(10) 


qQ QOi 1 qO QO 1 qO QO 1 qO QO 

•^jk'^lm "T •-^Im'^jk '^Ik'^jm '^jm'^lk 


jr Qa 1 e Qa , r Qa , jr Qa 
— '-'jk'Jlm + OlmOjk -t- OlkOjj^ -\- OjmOlk 


, (11) 


QO c/3 j_ qP QO I Qa c/3 _|_ C/3 Qo 
'-'jk'-'lm "t" '-'jk'-'lm '-'jm'-'lk '-'jm'-'lk 




(12) 


•^00^ Ql qO QO 1 qO Qa 
15 'Jjfe'Jim '-'jm'-'lk "r '-'ik'-'jm 




(13) 



These relations generate a long list of succeedingly less general relations by systematically 
equating all possible combinations of site indices. Setting j = k and / = m in eq. @, 
for example, implies that the operators ^S"" are spin density operators, 

[Sf,S^,] = 5,k2ie-^^S] . (14) 

The Hamiltonian (|ip now assumes the following form 

H = T.t^kS% + uY.{is^-i?-l) ■ (15) 

j,k 3 



6 



Since the particle number /° = J2j Sj is conserved, only the term (5*°)^ is relevant in the 
interaction part of the Hamiltonian. The other terms can be removed by a shift of the 
chemical potential. We retained them here to make obvious the invariance of H under 
the transformation The operators of the total spin are J" := ^J^jS"". It follows 
from that they commute with the Hamiltonian. Their s/2 commutation relations are 
obtained by summing (p!4D over j and k, 

[r,/^] =ie"^^r . (16) 

Now everything is prepared to formulate our main result. Consider the Hamiltonian 
(p!5|) with yet unspecified antisymmetric hopping matrix, tjk = —tkj- Let 

J"--=-2T.{ifj'^ + hAS^ + S'k-2))Sf, + g,ke''f'^S^S2) , (17) 

where gjk and hj^ are odd functions, and fjj = gjj = hjj = by convention. Then 
H commutes with J" if and only if the following functional equations between the 
coefficients are satisfied. 



{9ji - 9ki)hjk = ^hjihki , j ^k^ly^ j , (18) 

iUfjk/2ho + Qjkhjk = -\'^hjihki , j^k , (19) 

I 

Y^Ujihu- fkih.i) = , (20) 

tjk = h^hjk ■ (21) 



Here Hq is a free parameter which fixes the scale for J". The only solutions to these 
equations correspond to the cases of trigonometric and hyperbolic hopping amplitudes 
(0) under consideration. In the trigonometric case we find 

/,fc = , (7iA: = |ctg(7r(j-A;)/iV) , h^^ = i^i^'\'n{3 - k) /N) , (22) 

whereas in the hyperbolic case 

f^^ = -^^n-^) ' 9,k = \ot\i{K{j-k)) , = is\i-\K{j - k)) . (23) 
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ho has to be real in order for J" to be selfadjoint. We choose ho = —sm{TT/N) in the 
trigonometric case and ho = — sh(K) in the hyperbohc case. It is an unexpected fact 
that J" does not depend on U in the trigonometric case. Hopping part and interaction 
part of the Hamiltonian commute separately with J". 

One easily checks that our conserved operator J"" turns into known generators of 
Yangians in various limiting cases. In the nearest neighbour Hubbard limit of the hy- 
perbolic model (k oo) we obtain fj^ 6\j^k\,i/U, gjk sign(j — k)/2, hj^ — > 0. 
After a canonical transformation c^o- ^■^Cja we recover the Yangian generator of Uglov 
and Korepin 0]. In the limit f/ — > oo at less than half filling the model reduces to the 
"t-0" chain |]12[ with all states with double occupancies projected out from the Hilbert 



space. Because at half filling hopping is not allowed anymore in this limit, one can 
set Sj = 1, and recover the Yangian generator of the Haldane-Shastry chain |^ or its 
hyperbolic counterpart to leading order in (t/U)'^. 

Indeed J° itself generates a representation of a Y(s/2) Yangian. This is our second 
result. The spin operators J" and the conserved quantities J" satisfy the relations 

[I\j^] = cx,ur , (24) 

[[J\ J% [F, r]] + [[F, r], [I\ J^]] = -46{a^,,o.p,c,^u + ap„,^p,cx,u){r, I", (25) 

where 5 = — 1 in the trigonometric case, 5 = 1 in the hyperbolic case, and the further 
abbreviations 

CA^, := ie^^" (26) 
{a;i,a;2,a;3} ■= I ^i^j^k (28) 



have been used. Eqs. ( plf ) and (^) together with the defining relation (|T6D of the s/2 
Lie algebra are Drinfel'd's definition of the Y(s/2) Yangian Equation (|24D says that 
the J" transform like a vector representation of s/2, and is easily confirmed for our 
J". Since both equations (p^) and (p5|) are homogeneous, we could have introduced 



a deformation parameter /i^ on the right hand side of the Yangian Serre relation (25). 
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Since this parameter merely fixes the scale of J" and has no deeper physical meaning, 
we suppressed it here. We have confirmed (0) by direct calculation. The calculation is 
lengthy. Before we comment on it we formulate our third result. 

Under the transformation the generators transform into an independent 

set of generators /'°, J'^ of another representation of the Y(s/2) Yangian. The two repre- 
sentations commute, hence can be combined to a Y(s/2)©Y(s/2) double Yangian. Their 
commutativity is non trivial and is of dynamical origin, i.e. it relies on the functional 
equations (|I8|) - (^) between the coefficients that define J"", J'^ . 

To check the Yangian Serre relation, the original formulation ( pSf ) is rather inappro- 
priate. We used the following simplification in the s/2 case instead. Let 

K'^ ■= J7] _ 45(j/3)2j" . (29) 

Then a short but slightly tricky calculation shows that (^) is equivalent to the equation 

[J", K^] + [J^, = . (30) 

The left-hand side of (30) has a property that turns out to be very useful in practical 
calculations. It is traceless. Assume we are given an operator J", and we do already know 
that it transforms as a vector representation of sl2- Then this knowledge assures the 
identity [J", K'^] = 0. It is therefore sufficient to show that the left-hand side of equation 
(^) is proportional to 6"'^. This is a severe simplification, since the symmetrisation of 
the commutator produces a lot of terms proportional to S"'^, which can be neglected 
according to the above argument. The explicit expression for K"' in our case is 

j,k,i (31) 
+2A,,ig,i + ^?H)e"^^^^^7 + ih,ie-^^A,,S^, - A,,S^,^)iS} - 5,})} , 

where Aj^ := fjk + hjk{Sj +S^ — 2). To verify (|30|) , we used the following relations among 
the coefficients fjk, gjk, hjk in addition to their defining functional equations above. 

fjkigji - gki) = ^ifjihki - fkihji) , j ^ k^l^ j , (32) 

gjkgji + gkigkj + gijgik = , jj^k^ly^j , (33) 

4^/|. + ^'fc = S , j^k . (34) 
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The homogeneity of the lattices has not been used in the verification of the Yangian 
Serre relation in the bulk. However, it is necessary to guarantee the commutativity of 
J° with the Hamiltonian. This situation is similar to the case of the Yangian symmetric 
spin chains. Therefore we conjecture the existence of a Yangian symmetric long range 
Hubbard Hamiltonian on an inhomogeneous lattice. In analogy to the spin chain case 
||T3| the generator of its Yangian symmetry might be constructed by adding "potential 
terms" to the second order Yangian generator K"", eq. (pQ]). 

At this point we would like to emphasize that Yangian symmetry does not imply 
integrability. Nevertheless, we strongly believe that the models considered here are 
integrable, and are special cases of a more general integrable non Yangian symmetric 
model with elliptic hopping amplitudes. The proof of integrability would provide the 
basis for an understanding of the Haldane-Shastry chain and the nearest neighbour 
Hubbard model on a common ground. At the present state of knowledge these models 
appear rather unlike. The integrability of the Haldane-Shastry chain has been shown by 
exploiting a mapping to a related dynamical model |T^, whereas the integrability of the 
nearest neighbour Hubbard model follows from its connection to an integrable system 



of two coupled six- vertex models ||15| . We expect that a proof of the integrability of the 
non-nearest-neighbour Hubbard models will reveal a more generic structure. 

A first application of our new Y(s/2) Yangian will be the classification of the Jastrow- 
like eigenf unctions of the "t-0" model ||12|. For the system at finite on-site energy the 
situation is more complicated. Not even the ground state wave function is known. There 
is evidence that the wave functions are neither of Jastrow-type as for the Haldane-Shastry 
chain nor of Bethe Ansatz form as in case of the nearest neighbour Hubbard model flOl 
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